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Abstract 

NMR is the technique of election to probe the local properties of materials. Herein we present the results of density functional 
theory (DFT) ab initio calculations of the NMR parameters for fluorapatite (FAp), a calcium orthophosphate mineral belonging 
to the apatite family, by using the GIPAW method [Pickard and Mauri, 2001]. Understanding the local effects of pressure on 
apatites is particularly relevant because of their important role in many solid state and biomedical applications. Apatites are open 
structures, which can undergo complex anisotropic deformations, and the response of NMR can elucidate the microscopic changes 
induced by an applied pressure. The computed NMR parameters proved to be in good agreement with the available experimental 
data. The structural evaluation of the material behavior under hydrostatic pressure (from -5 to +100 kbar) indicated a shrinkage of 
the diameter of the apatitic channel, and a strong correlation between NMR shielding and pressure, proving the sensitivity of this 
i technique to even small changes in the chemical environment around the nuclei. This theoretical approach allows the exploration of 
all the different nuclei composing the material, thus providing a very useful guidance in the interpretation of experimental results, 
particularly valuable for the more challenging nuclei such as 43 Ca and 17 0. 
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1. Introduction 

Calcium orthophosphates, and in particular apatites (gen- 
eral formula Ca5(PC»4)3X (X=F,Cl,OH)), are extensively stud- 
ied because of their significance in many fields, from geol- 
ogy and materials science to medicine. Apatites, with vari- 
ous degrees of substitutions and defects, are commonly found 
in the Earth's crust and are considered one of the most abun- 
dant sources of phosphorous in the marine environment, act- 
ing as phosphorous "sinks" JT]. Apatites are also attracting a 
great interest in the field of biomaterials due to their intrinsic 
biocompatibility and bioactivity. Calcium phosphates, such as 
synthetic hydroxyapatite (HAp), fluorapatite (FAp), a and f3- 
tricalcium phosphate (a-TCP, /J-TCP), displayed intrinsic os- 
teoinductive properties, namely the capacity of promoting bone 
formation without presence of oesteogenic factors, as well as 
osteoconductive properties (the capacity of supporting and guid- 
ing the growth of newly forming bone). Apatite bioceramics are 
used in bone repair and augmentation, as bone cements, as coat- 
ings for bio-inert prostheses, and for bone scaffolding fl2j|3]. 
The study of natural apatites is complicated by the fact that the 
apatite crystal structure can readily accept substitution(s) from 
a large variety of ions. For example, fluoride ions are found 
to readily substitute in human enamel, due to their presence in 
drinking water as a result of soil erosion or added on purpose 
for caries prevention policies [4]. This capability of accepting 
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various substitutions and incorporating even bulky ions, such as 
Cd(II) and Cu(II), as well as U(IV), is making this class of ma- 
terials also attractive for nuclear waste management and water 
remediation [6). 

Fluorapatite (FAp) is one of the most well-characterized 
members of the apatite family, and in materials science it has 
found interesting applications as a host for solid state diode 
lasers when doped with rare-earth elements, such as praseodymium 
or ytterbium Q. It can be also used as a phosphor for fluores- 
cent lamps or plasma display panels, exploiting its photolumi- 
nescent properties when doped with Sr, or rare earth elements 

Good knowledge of the local structural properties and the 
response of FAp to external stimuli is of paramount importance 
to gain understanding of this material and, due to its structural 
similarity, the study of FAp can be considered as a reference 
for the other members of the apatite family. Solubility, crys- 
tallization, dissolution, and ion exchange have been intensively 
studied IfTOlfTTI . but relatively few studies explored the behav- 
ior of FAp under pressure and they were mainly focused on 
the bulk properties of the material, studied by using vibrational 
spectroscopic techniques and X-ray diffraction lfT2T - [T5l . 

In the present work we aim to explore the effects of pressure 
at the local level, by first principles calculations of the NMR pa- 
rameters response. NMR is the technique of election for study- 
ing the local structure of materials since it is very sensitive to 
changes in the atomic environment. A deformation occurring 
around a nucleus will in fact affect its shielding, thus providing 
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information on the deformation itself. In particular, our theoret- 
ical approach allows us to obtain unambiguously all the compo- 
nents of the shielding tensor and therefore to have a clear three 
dimensional picture of the effects of pressure on each nucleus 
in the material. 

The ab initio modeling of apatites is demanding from a 
computational point of view and was performed only in few 
cases in the past lfT6l - (2"7l . The published works have been mainly 
focused on structural and electronic features: a first principles 
study of fluorapatite and hydroxyapatite has been recently pub- 
lished by Rulis et al. (Ref. [25 1) where the authors computed 
electronic structure, charge distribution and X-ray absorption 
spectra on a supercell slab along the (001) surface. Some theo- 
retical NMR studies are available in the literature, but they are 
mainly focused on hydroxyapatite [28 - 30 1 . First principles cal- 
culations of NMR parameters have been available for molecules 
and clusters since the '70s ||3lT[35l while their application to pe- 
riodic/crystalline systems was implemented at a later time [ 36 1 
due to the difficulty of including a macroscopic magnetic field, 
requiring a non-periodic vector potential. In 2001, Pickard and 
Mauri implemented the Gauge Including Projector Augmented 
Wave (GIPAW) method [37] in the plane- wave pseudopotential 
(PWPP) framework, based on linear-response, aimed at calcu- 
lations of periodic systems. In the GIPAW method, the key to 
an accurate evaluation of NMR chemical shifts is the recon- 
struction of the all-electron wavefunction from a pseudopoten- 
tial calculation, by a modified PAW transformation, in order to 
describe the nodal structure of the wavefunction close to the 
nuclei. Hence the GIPAW method retains the accuracy of all- 
electron calculations at the much lower cost of pseudopotential 
calculations, i.e. explicit treatment of valence electrons. 

In the present work, the structural properties of FAp will 
be evaluated when hydrostatic pressure is applied in the range 
from -5 to +100 kbar, and particular focus will be dedicated to 
the response of the NMR shielding. For the nuclei with spin / > 
!/2 the quadrupolar coupling constant (C 9 ) and the asymmetry 
parameter {rj q ) will also be provided. This paper is organized 
as follows: in sections[2]and[3]we discuss the methodology and 
the structural data, in sections [4] and [5] we present and discuss 
our NMR calculations results, and conclusions are in section|6] 

2. Computational methods 

The calculations of fluorapatite (FAp) have been performed 
by using ab initio density functional theory (DFT) 0811 as im- 
plemented in the Quantum-ESPRESSO package |39|, with the 
PBE exchange-correlation functional 11401 . Well converged plane 
wave basis sets were employed for all the presented calcula- 
tions, using norm-conserving GIPAW pseudopotentials with a 
kinetic energy cut-off chosen at E cut - 80 Ry. The conver- 
gence criteria for the total energy was set in all cases at 10~ 8 
Ry. The k-points sampling in the Brillouin zone (BZ) used a 
Monkhorst-Pack grid with a converged mesh of 2 x 2 x 2. 

The FAp structure used for comparing the computed NMR 
parameters with the experimental data was constructed using 
the experimental lattice parameters (Ref. |41|) and fully re- 
laxed internal degrees of freedoms. The structures used for 



the evaluation of the NMR response under an applied pres- 
sure were obtained by imposing a target pressure (from -5 to 
100+ 1 kbar), and relaxing not only the atomic positions, but 
also the lattice parameters (Wentzcovitch damped dynamic al- 
gorithm, Ref. [ 42 J ) . The only constraint imposed on these re- 
laxations was the conservation of symmetry, under the assump- 
tion that no phase transition would occur at small enough pres- 
sures. This variable cell relaxation approach well simulates hy- 
drostatic pressure for bulk systems. 

NMR shieldings critically depend on the electronic wave- 
functions in proximity of the nucleus. However, those are ex- 
plicitly neglected in the pseudopotential approach, where the 
electronic wavefunctions are described by smooth functions and 
thus their correct nodal structure is lost. In the present work, 
specially designed norm-conserving GIPAW pseudopotentials 
were usecQ which allowed the calculation of the NMR shield- 
ings. We also used a modified calcium pseudopotential, con- 
structed with a rigid shift of the 3d orbitals by +3.2 eV, for the 
calculation of Ca and O NMR shieldings. The choice of using 
this pseudopotential was based on the observation by Profeta 
et al. [43 1 that in calcium oxides and calcium aluminosilicates 
the Ca PBE functional tends to overestimate the degree of co- 
valency of the Ca-0 bonds, thus causing the shielding on both 
Ca and O to be too small. However, this correction on the Ca 
pseudopotential had a negative influence on the computed equi- 
librium volumes and the computed phonon frequencies. There- 
fore in the present work, the structural relaxations and the NMR 
shieldings of F and P were obtained by using unit cells com- 
puted with the unshifted Ca pseudopotential, while the use of 
the shifted pseudopotential was limited to the calculation of Ca 
and O NMR parameters, assuming that the underprediction of 
the NMR shielding was a problem limited to Ca-O interactions, 
and that this effect was not relevant when calcium was present 
in different atomic environments. 

The NMR absolute chemical shielding tensor (*cr ) is de- 
fined as B md = -cr (r) • B ext , where B ext is the externally ap- 
plied uniform magnetic field. B ext generates a local current at 
the nucleus, which in turn induces a local magnetic field B md 
so that the total magnetic field experienced by the nucleus is 
B = B ext + B ind . 

Several different conventions for notation are in current use: 
in the present work the "Haeberlen notation" is used [44]. We 
report the isotropic NMR shielding (<x, JO ) as the average of the 
principal components of the shielding tensor (*cr ), expressed in 
a suitable set of axes X, Y, and Z: 

1 

o~iso = j(o-xx + ct Y y + o~zz), 
where the three principal components are related as follows: 

kzz - o- Lm \ > \o- xx ~ cr iso \ > \o- YY ~ cr iso \. 



1 We used standard Martins-Troullier pseudization. For oxygen and fluorine, 
the local channel is p and cutoff radii are 1.3 and 1.4 atomic units, respectively. 
For phosphorous, the local channel is d and cutoff radius is 1.9. For calcium, 
the local potential is derived from the d channel, and cutoff radii are 1 .45, 2.00, 
1 .45 for the s, p, and d channel respectively. 
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To fully describe the local symmetry around the nucleus, 
the shielding anisotropy (£) is also reported as: 

C = 0~ZZ - 0~iso\ 

as well as the asymmetry parameter (rj) as: 

(Tyy - crxx 



Finally, for the nuclei with / > V 2 ( n O and 43 Ca), the 
quadrupolar interactions are described by their quadrupolar cou- 
pling constant (C q ) and the asymmetry parameter (i] a ), defined 
as follows: 

r _ eVzzQ 
Lq ~ h 
Vxx - Vyy 

where e is the elemental charge, h is the Planck constant, and Q 
is the nuclear quadrupole moment (2( 17 0)= -2.56 x 10' 30 m 2 , 
2( 43 Ca)= -4.90 x lO" 30 m 2 ; Ref. El). V xx , V vv , and V zz are 
the diagonal elements of the electric field gradient (EFG) tensor 
in the principal axis system (PAS), with | Vzz \^\ Vyy 1^1 Vxx I- 
Most of the experimental NMR data are reported normal- 
ized in ppm as chemical shifts (6j SO ), with respect to a reference 
compound, where the chemical shift of the sample is obtained 
as: 

Osample — (X1U ) 

Vref 

where v samp i e and v re f are the resonance frequencies of the sam- 
ple and the reference compound, respectively. 

The chemical shift and the chemical shielding are related by 
the following: 

^sample — @~ ref ®~ sample 

To readily compare the computed data with the corresponding 
experimental values, it is necessary to calculate their chemical 
shift by choosing an appropriate reference material. Experi- 
mentally, most of the reference compounds for solid state NMR 
are molecules dissolved in water or an organic solvent, which 
are not well suited for calculations with our methodology. Our 
approach for overcoming this obstacle was therefore to choose 
solid reference compounds (secondary reference), with known 
experimental chemical shifts with respect to the standard liq- 
uid reference compounds (primary reference) recommended by 
IUPAC B6l . In Table [T] the solid reference compounds used 
in the present study are reported, along with their experimen- 
tal chemical shifts (<J^f) and their computed isotropic shielding 
values (o~™J). When more than one value for the experimental 
chemical shift was available, a range of computed NMR shifts 
is reported. 

The computed NMR chemical shifts (5, so ), for comparison 
with experiments, described in the present work are calculated 
as follows: 

r cr , ere/ f ref s ref 

disc = S iso + 6 iso = (cr iso - cr iso ) + 5 iso 



where 6' is the computed chemical shift with respect to the 
secondary reference material, while (5, JO is the chemical shift 
with respect to the primary reference material, which by def- 
inition has a 6 = ppm. The cr e ^ is the computed isotropic 
shielding for the secondary reference compound, while <t 1jo is 
the isotropic shielding computed for the nucleus of interest. 

Table 1: Primary reference materials (5 = ppm) and secondary reference 
materials with their experimental isotropic chemical shift (6™J with respect 
to the primary references | 47 |), and their computed NMR isotropic shielding 
(cr. ). 

v I SO ' 



Nucleus 


Primary 
reference 
Siso = (ppm) 


Secondary reference 

~re f ref 
5- CT- 
iso ISO 

Material (ppm) (ppm) 


43 Ca 


CaCl 2(a?0 (1.0M) 


CaO w 


136 


1032.79 


19p 


CFC1 3(/ ) 


CaF 2(l ) 


-107 
-104.8 


219.90 


31p 


H 3 P0 4(fl9 .) 85% 


BP0 4 


-29.5 
-31.2 


317.82 


17 


D 2 (/) 


quartz 


38 + 2 


211.84 



The proposed methodology of referencing the computed NMR 
shielding values allows for a direct comparison with most of the 
reported experimental data. Moreover it allows the referencing 
of our calculated values to systems which are computationally 
more similar to the ones studied in the present work, avoiding 
the complications inherent in simulating molecules in liquids. 

3. Structural data 

Fluorapatite has a hexagonal unit cell with P63/111 space 
group (n. 176), containing two formula units (42 atoms) per 
unit cell. The structure used in the present work is based on 
the crystallographic data of Sudarsan et al. (Ref. BP ). The 
nonequivalent atoms after relaxation of their positions are pre- 
sented in Table [2] The relaxed structure is very close to the 
experimental: the largest difference was of -0.12 A on the y 
coordinate of Ca(2), while all the remaining atomic coordinates 
were within +0.02 A of the experimental values. 



Table 2: Fluorapatite unit cell parameters and nonequivalent atoms obtained 
after relaxation of their atomic positions. Space group Pd-ifm (n. 176). 



Lattice 
Parameters 


Atom 


Wyckoff position 


X 


y 


z 


a=b=9.363A 


Ca(l) 


/ 


1/3 


2/3 


0.0006 


c=6.878A 


Ca(2) 


h 


0.2417 


-0.0054 


1/4 


a=/?=90 o 


P 


h 


0.3978 


0.3686 


1/4 


7=120° 


O(l) 


h 


0.3254 


0.4858 


1/4 




0(2) 


h 


0.5905 


0.4680 


1/4 




0(3) 


i 


0.3400 


0.2546 


0.0684 




F 


CI 








1/4 



The fluorapatite unit cell, and four adjacent FAp unit cells 
seen along the c-axis are presented in Fig.[TjA) and[TjB), re- 
spectively. Two distinct Ca sites are present in the unit cell: 
Ca(l) and Ca(2). Four Ca(l) occupy two vertical columns par- 
allel to the c-axis, and are usually referred to as "columnar". 
The remaining six Ca(2) can be divided into two groups of three 
calcium atoms, disposed to form two triangles lying on two 
planes perpendicular to the c-axis. The two triangles formed 
by the Ca(2) are respectively rotated by 60° and they delimit 
the apatitic channel where the fluorine ions are located. The 
presence of this channel is a feature common to the minerals 
with apatitic structure. Three nonequivalent oxygens are also 
present, and they are arranged to form six slightly distorted 
tetrahedra around the six equivalent phosphorous atoms. In par- 
ticular, each tetrahedron contains one O(l) and 0(2), and two 
0(3). The oxygen atoms form a network connecting Ca(l) and 
Ca(2) atoms, which are nine-fold and six-fold coordinated with 
oxygens, respectively. 

4. Computed NMR and comparison with experiments 

The crystal structure, as described in Table [2] was used 
for the NMR parameters calculations, and the results are re- 
ported in Table [3] Our calculated NMR chemical shifts are in 
good agreement with the experimental values available, espe- 
cially for 31 P where it is within 1 ppm. For the calcium and 
fluorine values, they are qualitatively similar to the experimen- 
tal data, within *13 and 8 ppm, respectively. All the nuclei 
showed a large shielding anisotropy, greater than 25 ppm (ab- 
solute value), with the exception of phosphorous atoms where 
it was about -8 ppm, thus indicating that the shielding around 
the P has a lower anisotropy due to the higher local symmetry 
of the phosphate group. 

According to the available literature, the only experimen- 
tal work on 43 Ca for FAp analyzed a fluoridated hydroxyapatite 
with a not specified fluorine content (Ref. [50|, w.r.t. CaCh 
sat.). In that work the two distinct calcium sites were not re- 
solved, although some asymmetry in the lineshape was evident. 
In Table [5] we report this measured 43 Ca chemical shift as the 
peak maximum, adjusted to a CaC^ (1.0 M) reference solution 
(+8 ppm). Throughout this paper all the reported 43 Ca chemical 
shift values are adjusted to this specific reference, which gives 
a 43 CaO chemical shift at 136 ppm (about 6-8 ppm higher than 
when referenced to a CaCl2 saturated solution). This approach 
does not totally eliminate the discrepancies in the calcium ref- 
erencing, since saturated solutions prepared from anhydrous or 
different hydrates of CaC^ can cause the 43 Ca chemical shift to 
cover a range of about 7 ppm, between +8.5 and +15.5 ppm. 
For a detailed description of how the calcium chemical shift is 
affected by differently prepared CaCl2 solutions, consult Ger- 
vais et al. (Ref. ||29l ). The few additional works regarding 
43 Ca are for hydroxyapatite, placing the 43 Ca chemical shifts at 
4.5-10 ppm and 17.5-22 ppm for Ca(l) and Ca(2), respectively 
1 29, 52, 53 1 . Even considering the possible discrepancies in the 
referencing, and the unspecified composition of the fluorapatite 
mineral analyzed in Ref. [50], these data suggest that when fluo- 
rine is present in the apatitic channel it causes an increase in the 
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shielding on the Ca nuclei (more negative chemical shift value), 
a trend which is well reproduced by our calculations. More- 
over, our calculations also showed that the presence of fluorine 
caused a smaller difference in shielding values between the two 
distinct calcium sites, when compared to hydroxyapatite. 

To the best of our knowledge, no experimental work was 
done on 17 for FAp, even though some data are available for 
hydroxyapatite, with n O chemical shifts at 108 and 115 ppm 
(unassigned) ll54l . A direct comparison of our data with these 
published values can potentially lead to errors, even though the 
crystal structures of the two compounds are similar. Neverthe- 
less, we do not expect large differences in the 17 chemical 
shift within the apatite class of compounds, since the environ- 
ment immediately around the oxygen atoms is scarcely modi- 
fied by different ions filling the apatitic channel. Our computed 
values for n O are in fact within approximately 1 1 ppm of the 
experimental values for HAp. Only two types of oxygens were 
detected in the experiment, even if three crystallographically 
distinct oxygens exist in the apatitic structure. From our cal- 
culations we found that the O(l) and 0(2) chemical shifts have 
very similar values, indicating a similar chemical environment, 
whereas the chemical shift of 0(3) is clearly different, and it 
is separated from the others by about 9-10 ppm. This observa- 
tion is in agreement with that found experimentally by Wu et al. 
(Ref. [54]) where the two peaks were separated by 7 ppm. If we 
approximate that the n O chemical shifts in FAp and HAp are 
similar, we can assign the more downfield peak found in HAp 
(115 ppm) to O(l) and 0(2), and the more upheld signal (108 
ppm) to 0(3). 

The quadrupolar coupling constants calculated in the present 
work for the two crystallographically distinct Ca nuclei had 
similar values of 2.3 and 2.9 MHz, but opposite sign. The ex- 
perimentally reported quadrupolar coupling constants for cal- 
cium (in HAp) are of 2.6 MHz for both Ca(l) and Ca(2) |52] 
l53ll . The difference in sign of the C q between the two distinct 
calcium nuclei found in our calculations was not observed ex- 
perimentally (at least in HAp), but it was found in calculations 
performed by Gervais et al. Il29l on HAp. All the computed C q 
for oxygens displayed a quite large and negative value of about 
-5 MHz. If our comparison with HAp holds, we can say that the 
the absolute value of the oxygens C q is similar to the ones ex- 
perimentally determined by Wu et al. (4.0 and 4.1 MHz) l54l . 

Although our calculated chemical shifts are in good agree- 
ment with the experimentally available values, some of the noted 
discrepancies can arise mainly from two facts: (1) variability 
in the preparation/composition of the primary reference com- 
pound, and (2) the error in the NMR parameters calculations. 
The variability in the preparation of the reference compounds 
becomes of greater influence for more challenging, low abun- 
dance and low sensitivity nuclei, such as 43 Ca and n O, where 
also fewer experiments are available to compare with. This er- 
ror can affect the value of both the primary and secondary ref- 
erence chemical shifts, thus generating a larger difference from 
the computed chemical shifts for those nuclei. For more abun- 
dant nuclei, where referencing is less challenging and more re- 
producible, the discrepancy between calculated and experimen- 
tal values is reduced. 



5. Structural evaluation and NMR calculations under pres- 
sure 

The FAp structure as described in Table [2] was fully relaxed 
without constraining either its lattice parameters, or its symme- 
try. The relaxed structure so obtained recovered all the orig- 
inal symmetries, while the resulting equilibrium volume was 
slightly bigger than the values reported in Table [2] (a = +1.6%, 
c= +0.22%). The tendency of overestimating the volume is a 
well known feature of the PBE functional. 

A set of structures were subsequently obtained by imposing 
a target pressure ranging from -5 to +100 kbar, and the com- 
puted energies and volumes for all these structures were used to 
calculate the bulk modulus of the material (Ko), by fitting with 
the Murnaghan equation of state. The computed bulk modulus 
was Kq =94 + 2 GPa, in good agreement with the experimental 
results reported by Comodi et al. (Ref. |i2|) of 93 GPa. The 
whole set of calculations at finite pressure did not cause ma- 
jor changes in the atomic positions within the crystal structure, 
nevertheless, when pressure was applied it caused a contrac- 
tion of the a and b-axis which was more pronounced than the 
c-axis. This result is also in agreement with experimental mea- 
surements of FAp under hydrostatic pressure, obtained using in 
situ single crystal X-ray diffraction, in the 0-7 GPa pressure 
range Ifl2l . This suggests that the c-axis is stiffer than a and 
b, possibly due to the arrangements of Ca(l) and F atoms into 
vertical columns along the c-axis. This columnar arrangement 
can cause a stronger electrostatic repulsion between the atoms 
within the columns, whereas the atoms that are staggered, such 
as Ca(2) and O atoms, can more freely rearrange in response 
to the applied pressure. The computed structures were further 
analyzed by calculating their pair distribution functions (PDF), 
using the ISAACS software [55|. This allowed us to explore 
the distribution of the atoms within the unit cell, and how the 
distances between the atoms changed when pressure was ap- 
plied. Most of the changes in the atomic distances occurred be- 
tween Ca(2)-F and Ca-O (Ca(2)-0(3) and Ca(l)-0(2)), sug- 
gesting a shrinkage of the apatitic channel and a closer packing 
of the phosphates, respectively. The P-O and O-O distances 
were essentially not affected by pressure, thus indicating that 
the phosphate units are more rigid than their surroundings, and 
confirming the larger degree of covalency of the P-O bonds. 

The NMR isotropic shielding (o", JO ) was computed and the 
results at the different pressures are reported in Figure [2] The 
shielding on the calcium atoms, in both sites, and on the fluorine 
atoms was found to decrease linearly with increasing pressure 
(linear correlation coefficient R 2 > 0.994). The applied pres- 
sure caused a contraction of the Ca(2)-F bonds from 2.33 A to 
2.23 A (4.3%) (linear fl 2 = 0.991). Whereas the change in bond 
length between Ca(2) and 0(3), the next calcium nearest neigh- 
bors, was from 2.36 A to 2.27 A (3.8%) (linear R 2 = 0.988). 
A similar trend of decreasing shielding on Ca atoms with de- 
creasing bond length was reported experimentally for a series 
of calcium containing compounds where the Ca 2+ were sur- 
rounded by oxygen atoms [29, 50]|56]|57). These authors also 
found that the deshielding was quite linear within a class of 
compounds (i.e. silicates, carbonates, aluminates, etc.), with 
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a common slope of 228 ppm/A. No clear correlation could in- 
stead be found by the authors when the Ca-X (X = halogen) 
were considered ll56l . A linear relationship between Ca shield- 
ing and Ca-O bond length was also observed in our calcula- 
tions for FAp {R 2 > 0.995), although the calculated slopes 
were larger (307 ppm/A and 414 ppm/A for Ca(l) and Ca(2), 
respectively). In the present work a clear linear relationship 
(R 2 > 0.996) was also found between the Ca shielding and the 
Ca-F distance, with slopes of 260 ppm/A and 352 ppm/A for 
Ca(l) and Ca(2), respectively. Also, in our case, the larger slope 
for Ca(2) shielding vs. bond length is a further indication that 
the applied pressure is mainly affecting the diameter of apatite 
channel, which is delimited by the Ca(2) atoms. 

Another point to be noted is that the FAp structure with- 
out applied pressure had a higher shielding on Ca(l) than on 
Ca(2), but at the highest pressure computed (100 kbar) the two 
shielding values were switched. This observation is somewhat 
in agreement with the fact mentioned by Lin et al. (Ref. |56|) 
where they observed that when the electronegativity of the atom 
connected to the Ca increased, the shielding on the calcium de- 
creased. Similarly here, when the distance between calcium 
and fluorine reached below a certain value, the very electroneg- 
ative F caused the shielding on the Ca(2) to decrease faster than 
the one on Ca(l), for which the first nearest neighbor is instead 
oxygen. The value at which this "inversion" was observed is 
2.28 A and it is just slightly smaller than the sum of the Ca 2+ 
and F~ ionic radii (2.30 A). This also suggests that the ionic 
radii of O and F may play an important role in determining this 
observed behavior: oxygen has in fact a slightly larger ionic 
radius (1.38 A) than fluorine (1.30 A), so that the effects of flu- 
orine on the shielding of Ca(2) become more relevant only at a 
shorter distance 

Interestingly, while the shielding on the fluorine atoms fol- 
lowed a trend similar to the calcium atoms, phosphorous dis- 
played an opposite behavior, with its shielding increasing with 
increasing pressure. As mentioned above, the change in P-0 
bond length was minimal in the range of pressures evaluated, 
thus even only a 9 ppm change in the shielding was a large 
change with respect to the change in bond length. This suggests 
that the shielding response of P is more susceptible to changes 
in its surrounding environment, and that this larger suscepti- 
bility together with the direct relationship between shielding 
and pressure is the result of the more rigid/covalent nature of 
the P-O bonds, with respect to the Ca-F or Ca-O. Finally, the 
shielding around the oxygen atoms was not greatly affected by 
pressure. This indicates a counterbalancing of the two opposite 
tendencies of increasing and decreasing of the shielding, when 
applying pressure, due to the fact that oxygens are bound to 
phosphorous while also bridging the calcium atoms. 

The pressure had also an effect on the absolute value of the 
shielding anisotropics, with an increasing anisotropy around Ca 
and F nuclei, and decreasing values for O and P. The latter find- 
ing is in agreement with the gradually decreasing difference be- 
tween the P-O bond lengths of the phosphate tetrahedra, indi- 



2 The ionic radii used are from Ref. 1581 . considering a coordination number 
C.N. = 3 for F. C.N. = 6 for O, and C.N. = 6 for Ca(2). 



eating a more symmetric environment around the P atoms. A 
decreased distortion of the phosphates from tetrahedral symme- 
try, with increasing pressure, was also observed experimentally 
by noticing a reduced site group splitting of the Raman signals 
of the V3 and V4 phosphate vibrations |15|. The different con- 
tributions to the shielding were also evaluated and they were 
mainly due to paramagnetic currents tightly localized near the 
nuclei, except for the anions (F and O) where the induced cur- 
rent was more delocalized towards neighboring atoms, where 
partially-covalent bonds are formed. Finally, the quadrupolar 
parameters were only minimally affected by the pressure, with 
a change of less than 0.6 MHz for C q and less than 0. 1 1 for r\ q . 

6. Conclusions 

We successfully computed the NMR chemical shifts and 
quadrupolar parameters of fluorapatite (FAp) and they resulted 
in good agreement with the available experimental data. An 
assignment was suggested for the 43 Ca (FAp) and n O (HAp) 
experimental chemical shifts (Ref [50 1 and [54|, respectively). 
The behavior of FAp under pressure and its effects on the NMR 
parameters were also successfully explored. The pressure caused 
a shrinkage of the apatite channel, with the largest changes in 
bond lengths observed between Ca-F and Ca-O, whereas the 
P-O bonds were minimally affected. A linear behavior of the 
NMR shielding vs. pressure was found in all nuclei, although 
Ca and F had an opposite slope with respect to P, while the 
change in the shielding on the oxygens was minimal. We be- 
lieve that the opposite behavior of Ca and F when compared 
to P is due to the different covalency/ionicity character of the 
bonds involved, with the more covalent bonds having a larger 
shielding with decreasing atomic distance, and vice versa. Fi- 
nally, our calculations can guide the calibration and the inter- 
pretation of the chemical shift of apatites in NMR experiments 
under high pressure. 
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Figure 1: (A) Fluorapatite unit cell and (B) four adjacent fluorapatite unit cells 
depicted along the c-axis. The channel (shadowed region) and the columnar 
calciums are indicated. (Ca = blue, F = green, P = purple, O = red; the images 
were obtained using VESTA 1481) 
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Figure 2: The NMR isotropic shielding (cr, so ) for the different nuclei in fluora- 
patite (FAp), when hydrostatic pressure was applied from -5 to 100 kbar. 
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